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Important E2 contributions to the (e,e') cross sections of low-lying orbital Ml excitations are 
found in heavy deformed nuclei, arising from the small energy separation between the two exci- 
tations with I^K = 2^1 and respectively. They are studied microscopically in QRPA using 
DWBA. The accompanying E2 response is negligible at small momentum transfer q but contributes 
substantially to the cross sections measured at 9 ^ 165° for 0.6 < q^H < 0.9 fm"^ (40 < Ei < 70 
MeV) and leads to a very good agreement with experiment. The electric response is of longitudinal 
C2 type for 9 < 175° but becomes almost purely transverse E2 for larger backward angles. The 
transverse E2 response remains comparable with the Ml response for geff > 1-2 fm~^ {Ei > 100 
MeV) and even dominant for Ei > 200 MeV. This happens even at large backward angles > 175°, 
where the Ail dominance is limited to the lower q region. 
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I. MOTIVATION 

Inelastic electron scattering at backward angle 9 = 
180° was used by Barber, Peterson et al. more 
than 30 years ago to study nuclear magnetic dipole (Ml) 
excitations at the linear accelerator in Stanford. The 
predominance of Ml excitations in backward scatter- 
ing is expected from qualitative considerations ^ of the 
(e, e') cross section in the plane wave Born approximation 
(PWBA). The cross section can be decomposed in this 
case into longitudinal (Coulomb) and transverse (elec- 
tric and magnetic) terms, multiplied by corresponding 
kinematical factors. The electric field of the incoming 
electron can have both longitudinal and transverse com- 
ponents, while its magnetic field is purely transverse. If 
the electron rest mass and the nuclear excitation energy 
can be neglected in comparison with the incident electron 
energy, which is often the case with low- lying excitations, 
the longitudinal kinematical factor vanishes for — 180° 
and only transverse multipoles are excited by the inelas- 
tic scattering 

The (e, e') cross section can be related further to the 
transition probability for photoabsorption through ap- 
proximate considerations within the "virtual photon" 
method. One arrives in this way to the rough quali- 
tative estimate |^ for dominant magnetic over electric 
transverse excitations (both of the same multipolarity) . 
This property is due to the momentum transfer during 
electron scattering. It stands in contrast with the corre- 
sponding relationship for radiative excitation, where no 
momentum is transferred and the EL excitation is one 
order of magnitude stronger than the ML excitation {L 
stands for multipolarity). 

A more definite qualitative conclusion about the mag- 
netic dominance at backward angle can be drawn by 
assuming a purely spin- flip transition Q]. The trans- 
verse electric field of the incoming electron is negligibly 
small in this particular case if the excitation energy 
is small in comparison with the transferred momentum 
q, i.e. Ex <C qfic. The same result is obtained for a spin- 
flip transition in the limit of small momentum transfer, 
q ^ £'a;/?ic All these considerations are made in the 
long- wave limit, qR <ti 1, where R is the radius of the 
target nucleus. The momentum transfer has throughout 
this article the inverse dimension of length (fm^^), i.e., 
q/h is denoted by q. 

The peculiar selectivity of the backward (e, e') scat- 
tering towards magnetic excitations was used extensively 
in the last decades to study Ml transitions mainly in 
spherical nuclei |^,^. The Ml operator can excite there 
particle-hole states, which are spin-orbit partners. These 
excitations are therefore often called spin-flip, though the 
orbital part of the Ml operator has also non-vanishing 
matrix elements In addition to these predominantly 
spin transitions, a qualitatively new type of mainly or- 
bital Ml excitations with I^K = 1+1 become possi- 
ble in deformed nuclei [Q. They are closely related to 



the so-called scissors mode, predicted by the two-rotor 
model [HH]. It describes isovector rotational oscillations 
of neutrons against protons, which exist only in deformed 
nuclei. Their experimental study started a decade ago 
through backward inelastic electron scattering at the lin- 
ear accelerator in Darmstadt, see e.g. the review articles 

We have calculated within the quasiparticle random- 
phase approximation (QRPA) the (e, e') form factors of 
low-lying orbital Ml excitations in Titanium | [Tl[ |, rare- 
earth |]l^, and actinide nuclei. The present study 
was motivated firstly by noticing at higher momentum 
transfer some discrepancies between our theoretical Ml 
form factors and experimental data. These deviations are 
not present in lighter nuclei but seemed to occur system- 
atically in heavy nuclei. The first diffraction maximum 
of the experimental (e, e') cross section is usually well de- 
scribed theoretically when the experimental B{M1) value 
is reproduced in QRPA. In contrast, the g-dependence of 
the theoretical cross section becomes very sensitive to 
the details of the microscopic wave function at higher 
momentum transfer q. 

Our theoretical Ml transition density of the 1+ exci- 
tation at 3.78 MeV in ^^Ti ^ describes well the corre- 
sponding experimental form factor p4| not only at low 
but also up to the highest measured momentum transfer, 
9cff < 2 fm^^. We were even able to predict ||ll| the form 
factor of the strongest Ml excitation (at 7.2 MeV) in this 
nucleus in the same q-range before the corresponding ex- 
perimental data became available |^5|. The comparison 
with experiment can be seen in Fig. 3 of Ref. Q. 

On the other hand, our theoretical Ml form factors 
for rare-earth nuclei p^ , though lying within the ex- 
perimental error bars, seem to underestimate system- 
atically the experimental values even at lower momen- 
tum transfer (for q > 0.6 fm^^). In contrast to ^'^Ti, 
high-g experimental data are not available for single 
Ml excitations in heavy deformed nuclei. We were in- 
clined to attribute initially the above discrepancies to the 
adopted parametrization of the Woods-Saxon potential. 
The K'^ — excitations in deformed nuclei are weakly 
collective fl^Jl^ , because the main collectivity is concen- 
trated in the spurious isoscalar rotational state. Thus, 
the microscopic wave functions of the low-lying orbital 
Ml excitations are very sensitive to the details of the 
single-particle level scheme around the Fermi surface. In 
order to improve the agreement with the experimental 
form factors, we decided to shift in rare-earth nuclei the 
energy of several single-particle levels close to the nuclear 
surface ||l^] , a procedure used formerly in microscopic cal- 
culations for spherical nuclei. The modified level scheme 
did not produce, however, a substantial improvement. 

Thus, we turn our attention to possible contributions 
from other multipolarities, which should have the prop- 
erty of appearing only in heavy nuclei. We are going 
to present in the next section preliminary qualitative es- 
timates for the expectation of E2 contributions to the 
(e, e') cross sections at backward angle. The microscopic 



2 



calculation of transition densities within our QRPA for- 
malism is described in Sect. III. The theoretical (e, e') 
cross sections, obtained from these transition densities 
in the distorted wave Born approximation (DWBA) for 
9 — 165°, are compared to experiment in Sect. IV. A 
qualitative analysis of our DWBA results for larger back- 
ward angles is presented in Sect. V using more transpar- 
ent PWBA estimates. The conclusions are summarized 
in Sect. VI. The most important ingredients of our for- 
malism can be found in Appendices A-D. 



II. PRELIMINARY ESTIMATES 

Our search for possible contributions from higher mul- 
tipolarities was stimulated by the experimental study 
of the E2 transition accompanying the strongest M\ ex- 
citation, = 3.075 MeV, in ^^^Gd. The E2 transition 
excites the first member (with I'^K ~ 2+1) of the rota- 
tional band built on the top of the vibrational state with 
I'^K — 1+1. This is the only accompanying E2 tran- 
sition studied until now in heavy deformed nuclei. The 
energy separation between the two excitations {M 1 and 
E2) corresponds to the rotational energy 
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where J is the moment of inertia of the rotational band 
developed on the intrinsic 1+ state. Since its value 
is not known, one could use as a rough estimate the 
ground-state moment of inertia, determined from the 
energy E[2X) = ih^/Jg.s. of the first rotational 2+ 
state in ^^^GA fl^. The so obtained energy separation, 
Ai?i2 — 59.3 keV, should be considered as an upper hmit, 
because the rotational bands of excited states have usu- 
ally larger moments of inertia than the g.s. band [ pO[ . 
The 2+1 state found in ^^^Gd [|8) lies indeed only 21 
keV above its 1+1 band head with E^ = 3.075 MeV. 
This energy separation is comparable with the energy 
resolution l^E^xj}. of the considered (e, e') experiments, 
20 keV < AEexp. < 40 keV, depending on the incident 
electron energy pl[ |. 

The largest backward angle was 9 = 165° in these ex- 
periments. It corresponds to a ratio between the trans- 
verse and longitudinal kinematic factors [Q 
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Although the longitudinal Coulomb excitation is strongly 
damped by this large factor, it could still provide a C2 
contribution to the cross section comparable with the Ml 
contribution. This is due to the fact that the longitudi- 
nal E2 transition matrix element, before being reduced 
by the above factor, is usually much larger than the Ml 
matrix element. A new 180° electron scattering facility 
is already in operation in Darmstadt p^. Though the 



aperture of such facilities leads to a somewhat smaller 
effective mean angle of 178.5° the longitudinal con- 
tributions are already negligible at this angle, where the 
ratio (H) amounts to 5848. On the other hand, the wors- 
ened energy resolution, typically AE^^p. — 100 keV at a 
fully backward angle ||23| , makes quite plausible the ap- 
pearance of transverse E2 contributions to the measured 
Ml cross section, if both are of comparable magnitude. 

The ground state of even-even nuclei corresponds to 
I'^K = 0+0, so that due to the selection rules each one 
of the above considered 1+1 and 2+1 states can be excited 
only by Ml and E2 transitions, respectively, but not by 
both of them. However, the multipolarity of the excited 
states studied is usually not known in experiment. Thus, 
most of the former qualitative estimates for the relative 
CL, EL, and ML contributions to the PWBA cross sec- 
tion consider transitions with the same multipolarity L 
in order to find out what kinematics would enhance se- 
lectively the ML excitation. 

In contrast, we are interested in the present study in 
the relative E2 and Ml contributions. They excite two 
different states, but contribute nevertheless to the exper- 
imentally measured cross section only because in heavy 
deformed nuclei the energy separation of the two states 
is comparable with the experimental energy resolution. 
This situation does not occur in lighter even-even nuclei 
whose small moment of inertia produces a large energy 
separation between the two states and the two cross sec- 
tions are easily resolved in experiment. We denote by 
E2 both longitudinal C2 and transverse E2 components, 
unless one of them is specified explicitly. 

Our study was motivated also by one of the first cal- 
culations of (e, e') cross sections in heavy deformed nu- 
clei iQ. The PWBA cross sections for members of the 
ground-state and /3- vibrational bands in ^^^Sm are ob- 
tained in the collective Tassie model assuming an incom- 
pressible and irrotational fluid. Although a purely longi- 
tudinal (Coulomb) excitation is described by this simple 
model, the qualitative g-behaviour of the calculated cross 
sections supports our expectations of improving the Ml 
cross sections by taking E2 contributions into account: 
i) the diffraction maxima of higher multipolarities have 
smaller amplitudes but are broader, ii) the I^K = 2+0 
state of the /3-vibrational band has even a larger cross 
section than the band head (0+0), i.e. the magnitude or- 
dering is inversed for the intrinsic excitation with respect 
to the g.s. band. These qualitative features were con- 
firmed experimentally some 20 years later |^ through 
(e,e') on ^^^Gd. 

The above arguments about possible appreciable E2 
electroexcitations even at backward angles stimulated us 
to examine more carefully the qualitative estimates, on 
which the commonly expected Ml dominance in back- 
ward electron scattering is based. The relative contribu- 
tion of the transverse magnetic and electric terms to the 
(e, e') cross section is given by the approximate PWBA 
relationship ||^: 
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where we have used the correct ratio from Eq. (30) of 
Ref. [^. The factor 10 is a rough estimate for the ratio 
between the two cross sections at the photon point where 
they are proportional to the respective reduced probabil- 
ities for radiative transitions. Thus, one should expect 
predominantly magnetic excitations in backward electron 
scattering, because in most experiments E^: <C qhc. But 
the qualitative relationship was derived within the 
"virtual photon" method valid for q — > Ex /he. It was 
stated explicitly |^ that the relationship (||) is not enough 
to upset in general the dominant magnetic character of 
backward electron scattering. More detailed considera- 
tions will be discussed in Sect. V below. 

Without recurring to virtual photons, Willey arrives to 
a more general qualitative relationship between the two 
cross sections l26|: 
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The magnetic transition is obviously dominant, because 
is derived in the long- wave limit, qR <C 1. However, 
the long-wave approximation is not fulfilled in the above 
discussed experiments with heavy nuclei [p^ , ^ , where 
q > 0.2 fm-i, R « 6.5 fm, and qR > 1.3. With this 
numerical value one could conclude from (|^) that the 
transverse electric should be comparable with the trans- 
verse magnetic electroexcitation, and even dominant for 
q > 0.5 fm~^, i.e. in the region where we expect impor- 
tant E2 contributions. Of course, (Q) is no more valid for 
such large values of qR. 

The above simple considerations are not able to pro- 
vide a reliable prediction for the relative interplay of elec- 
tric and magnetic electroexcitations in backward scatter- 
ing. They ignore the dynamics contained in the wave 
functions of the studied nuclear excitations. More defi- 
nite conclusions can be drawn only from realistic micro- 
scopic calculations which are able to provide satisfactory 
quantitative description of the experimental data. We 
are using the DWBA because the large charge of heavy 
nuclei produces a strong Coulomb potential distorting 
the electromagnetic field of the incoming electron. After 
our first encouraging results on ^^^Gd p7p8], we are go- 
ing to study here the problem of electric contributions to 
backward scattering more systematically and to compare 
our theoretical results with all the available experimen- 
tal cross sections for single Ml excitations in rare-earth 
nuclei. 



III. TRANSITION DENSITIES IN QRPA 

We solve the QRPA equations of motion for intrinsic 
excitations with K'^ — I"*' using a model hamiltonian 



F{m,t) = [Ho,J{m,t)], t ^ n, p, r = (5) 
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where the isospin index t — n,p denotes neutrons or pro- 
tons. The quasiparticle (q.p.) mean field Hq ^6 30 
is given by an axially-symmetric Woods-Saxon poten- 
tial H plus BCS pairing. Apart from the deformation, 
all the remaining Woods-Saxon parameters are given for 
each nucleus by the isospin dependent parametrization 
of Ref. [Hj. It is based on a single set of isospin inde- 
pendent parameters for all nuclei without respect to their 
mass and charge numbers. The hexadecapole deforma- 
tion is taken from experiment psf , while the quadrupole 
one is determined for each nucleus by fitting the micro- 
scopically calculated ground state quadrupole moment 

|. As in 
ator con- 



to the corresponding experimental value |^ 
most of our previous works, we use an oscil 
stant hcjn — 50 A^^^^ MeV and include in the basis (Al) 
all the os cilla tor eigenfunctions with principal quantum 
numbers ( |A^ ) J\f < 11. The pairing gaps are chosen to 
reproduce roughly the even-odd mass differences P5[ . 

The next two terms in the model hamiltonian ^ are 
separable residual interactions of quadrupole and spin 
type, respectively. The former interaction is constructed 




derived from 
id^. The to- 



with operators of quadrupole type F 
the deformed part of the q.p. mean field 
tal angular momentum J(m, t) ([ij) is symmetrized over 
the signature m = ±1, see Appendix B. The isoscalar 
coupling constant kg of -ffpF is determined microscop- 
ically |3q| in order to restore the rotational invariance 
of the model hamiltonian, violated by the deformation. 
The value r = — 2 is adopted for the ratio between 
the isovector and the isoscalar coupling constants. This 
choice, discussed in |29 , gives rise in QRPA to a high- 



lying E2 strength, whose energy distribution reproduces 
qualitatively the isovector giant quadrupole resonance 
observed in (e, e') experiments. We adopt the above 
rough round value, because of the uncertainties of the 
scarce experimental data on heavy deformed nuclei. Our 
attention is focused here on the low-lying orbital M 1 exci- 
tations, which are not very sensitive towards the strength 
of the isovector quadrupole interaction. They are in- 
fluenced much stronger by the exclusion of spuriousity 
through restoration of the rotational invariance. They 
acquire afterwards a well pronounced character of isovec- 
tor rotational vibrations [p6 17|. 

The last term in is the spin-spin interaction Hss 
prt. Its coupling constants. 



c(+)A = 200 MeV, c(-) = -0.5c(+), 



(6) 



are derived from nuclear matter calculations using the 
Reid soft core interaction in the way explained in [ p^ . 
The coupling strengths (H) provide a good description 
p?! of the experimental spin Ml strength in rare-earth 
nuclei, deduced from inelastic proton scattering. 

We present here results for all the rare-earth nuclei 
with experimentally measured (e, e') cross sections of sin- 
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gle 1+ excitations. They are characterized by a relatively 
strong Ml transition and an excitation energy of about 3 
MeV. The energies E.^; of the Ml excitations, the B{M1) 
and B{E2) values, and the orbit-to-spin ratio Ro.s. of 
the Ml matrix elements are compared with experiment 
in Table ||. Expressions for Ro.s., Ml and E2 transition 
probabilities in terms of the RPA amplitudes are given 
in The B{M1) values for transitions to I^'K = 1+1 
states are calculated with bare orbital and effective spin 
gyromagnetic factors, gs = 0.7gl^°°. The B{E2) values 
for transitions to the corresponding I'^ K = 2+1 states, 
discussed already in Sect. II, are obtained without intro- 
ducing effective charges, i. e. they are of purely proton 
nature. The same orbital and spin gyromagnetic factors 
and effective charges are used, of course, in the calcula- 
tion of transition densities, presented below. 

It is seen from Table | that the theoretical results 
are in a good agreement with the available experimen- 
tal data from (e,e'), ip,p'), and (7,7') scattering. One 
should note especially the agreement with the experimen- 
tal B{E2) value for ^^^Gd, extracted from (e,e') [|l|. It 
is closely connected with the main topic of the present 
study: the E2 response accompanying the Ml electroex- 
citations in backward scattering. Moreover, as discussed 
in Sect. II, this is the only accompanying E2 transition 
identified experimentally as yet. 

The experimental situation is more complicated in 
issQd, where the (e,e') peak || at 3.2 MeV with 
B{M1) = 1.4(3) is resolved in nuclear resonance fiuo- 
rescence [|l] in two Ml transitions: the first one is at 3.19 
MeV with B{M1) = 0.66(8) /^^, and the second one is 
listed in Table |. As discussed in Sect. II, each Ml tran- 
sition is accompanied by a corresponding E2 excitation. 
The latter was studied experimentally only in ^^^Gd as 
yet, where it has an excitation energy £^2,(2+1) — 3.096 
MeV |l|. 

The large ratio -Ro.s. 0] between the orbital and spin 
Ml matrix elements of each excitation, listed in the last 
column of Table |, indicates that these are predominantly 
orbital transitions. Most of them represent the strongest 
single Ml transition in the corresponding nucleus. They 
were analysed by comparison with the microscopic re- 
alization [IlSi of the scissors mode H, as well as with 



12|. This analysis allows 
y collective scissors mode 



synthetic rotational phonons [ 
us to interpret them as a weak 

Electroexcitations are described by the following one- 
body charge and current density operators [^,||: 



P(r) 



A 

• c 
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J(r)=r(r)+f(r), 



j^(r) 
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(7) 

(8) 
(9) 



A(r) 



2Mc 
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where e,M, e and g'^ are the unit charge, the nucleon 
mass, effective charges and spin gyromagnetic ratios, re- 
spectively. p{y) and J(r) are scalar charge density and 
vector current density operators, respectively. The latter 
is a sum of the convection current j^(r) and the mag- 
netization (or spin) current j^(r). The summation runs 
over the coordinates of all the nucleons. 

The momentum operator p in (^ has to be hermitian 
symmetrized in order to satisfy the continuity equation 



sym. j_ 



(P + P^)- However, we do not need this 



symmetrization in (||) because we work always with sig- 
nature symmetrized operators, defined below (p^. They 
possess even higher symmetries, listed in Table |l|. Ap- 
pendix B, including the hermitian property required for 
the convection current (^. 

In consistence with our QRPA calculations, we do not 
use effective charges in (0,||), but the spin operator is 
renormalized by a factor of 0.7, included in the spin gy- 
romagnetic ratios: 

£" = 0, e^' = 1, 5,' = 0.7g,^(free), 
g,^(free) = -3.8263, <7p(free) = 5.5858. (11) 

Thus, only protons contribute to the calculated transition 
densities of charge and convection currents. 

The charge and current density operators (|7[||) can 
be expanded in spherical harmonics and vector spheri- 
cal functions, respectively Q]: 



LM 

J(r) = ^ Yl^,,,{mLL'M{rl 

LL'AI 

YLL'Aim = J2(L'fi, lfi'\LM)YL'^in)e^,. (12) 
mm' 

The resulting multipole density operators have the form 



PLMir) = J Ylm{^)p{^) dil, 

JLL'Aiir) = J YLL'Min)>J{r) dn, 

L' ^L.L±1. 



(13) 



f{r) = V X A(r), 



We are interested only in Ml and E2 electroexcita- 
tions. The former are generated by the transverse Ml 
current density operator J7ii'm('')- The latter receive lon- 
gitudinal (Coulomb) contributions from the charge den- 
sity operator p2M{r) and transverse contributions from 
the E2 current density operators J'2L'm{t), L' = 1,3. 
Only operators with M = ±1 contribute to the consid- 
ered intrinsic excitations with = 1+. 

We use m-symmetrized operators in our QRPA for- 
malism [ff6|. This is necessary to close the algebra of 
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the commutation relations. Such operators possess also 
definite symmetries with respect to complex, hermitian, 
and time conjugation by simply changing their phase un- 
der these transformations, see Table || in Appendix B. 
We are able to ensure in this way the hermiticity of the 
convection current (^) , to take further symmetries of the 
microscopic expressions analytically into account, and to 
reach the minimal set of irreducible terms. The signature 
m = ±1 corresponds to the intrinsic x- and y-axes, which 
are indistinguishable in axially symmetric nuclei. Thus, 
final expression for physical observables are independent 
of m. 

There are two main types of operators in the quasibo- 
son formalism: longitudinal and transverse, according on 
their phase a — ±1 (B7) with respect to combined her- 
mitian and time conjugation, considered in Appendix B. 
Operators of the same type co mm ute with each other in 
the quasiboson approximation (B9). Non-commuting op- 
erators are, e.g., the quadrupole operator Q{m) and the 
total angular momentum J(m) being of longitudinal 
and transverse type, respectively. We define here addi- 
tionally m-symmetrized multipole charge p2{'m,r) and 
current jLL'ifn^r) density operators using ([l3|): 

J(m) = i[mJ+ - J±^X±iJy, 

Q{m) = \[mQ2i + Q21]: 

P2{m,r) = i[m/52i(r) 

jLL'{m,r) = ^[mJLL'iir) + jlLai^)]. (14) 

It is seen from their quasiboson representation (B5) that 
the charge density is of longitudinal (quadrupole) type 
with a — +1, while the current densities are of transverse 
(momentum) type with cr = — 1. The definitions jl^ ) 
ensure that the quasiparticle (q.p.) matrix eleme nts o f 
all transition operators have the same symmetries (Bll). 

The (e, e') transition densities are the reduced matrix 
elements (r.m.e.) of the multipole density operators ( [T^ ) 
between the initial and final nuclear states We de- 
scribe here in QRPA intrinsic excitations with K"^ = 1+ 
in axially symmetric even-even nuclei, where K is the 
projection of the total angular momentum /, on the in- 
trinsic symmetry axis of the nucleus and tt is the parity 
of the wave function. The states are characterized by 
I'^K, ly, where = 1, 2, ... labels the intrinsic QRPA ex- 
citations rj,(m)|) (pl|), obtained by solving the QRPA 
equations of motions. 

The transition TL with multipolarity L (Ml or E2) 
takes place between the ground state with I'^K — 0+0 
and an excited state ly with I^K = L+1. The transition 



densities (C7), correspo ndin g to the multipole operators 
(ISl), are obtained from (C4). We use the final expression 
(C9) in terms of q.p. matrix elements (CS}of step-up 
operators, because the microscopic sum (C9) contains 
the minimal set of independent terms, not related with 
each other through any symmetry. Thus, the E2 and Ml 
transition densities for the v-th QRPA excitation have 
the form: 



p^(r) = (2n,Hl4(^)llO) = i 

0<i<f 

Fl^[fi)p2i{fi, r) + F;i(/z)p2i(/1 r)l , (15a) 



o<i<S 



JI^^,ir)^z{Ln,,y\\jl^,ir) 



F'^,if^)JLL'lif^,r) + F'^,{f^)JLL'liflr) 

JEL'ir) = JLL'ir) + JL'L'ir), (15b) 

where p2i(.f», r), JLL'i ifi, r), ar e the q.p. matrix ele- 
ments (D3), (|D4|), (D7), and (D8) of the transition den- 
sity operators P21M (longitudinal) and Jhuiir) (trans- 
verse), given by (|l3|). The transverse transition density 



( |l5q ) is decomposed into a sum of convection and spin 
(or magnetization) parts, according to (|8|)-([l0|). The in- 
formation about the considered i^-th QRPA excitation is 
containe d in the factors F^{fi), which a rc li near combi- 
nations ( |C10| ) of the QRPA amplitudes ( p] ). 

The QRPA transitions densities for Ml and E2 exci- 
tations of the K"^ = 1+ state in ^^^Gd from Table | are 
displayed in Fig. |l|. This is the strongest low-lying (or- 
bital) Ml transition in ^^^Gd. The longitudinal charge 
transition density P2{t') has the largest amplitude (note 
the different ordinate scales), followed by the transverse 
Ml transition density Jii{r). Each transverse density 
(continuous curve) is a sum of three components, aris- 
ing from the proton convection (long dashed curve) and 
spin (dotted curve for protons and short dashed curve for 
neutrons) currents (||). In the case of Jii{r) and 1/21 (r) 
the proton convection current is dominant, while the spin 
currents for protons and neutrons are out-of-phase, as it 
should be expected for an isovector vibration, and cancel 
each other to a great extent. 

We do not display the transition densities of the other 
1+ excitations from Table ^, because they are very similar 
to those in Fig. ^ apart from the cases of ^^*Gd and 
^^^Er. Since the latter two cases exhibit also similarities, 
we show only the transition densities for ^^^Er in Fig. ||. 
This excitation has a more pronounced volume character 
because, in contrast to Fig. 0, the transition densities 
are no more peaked near the nuclear surface. Although 
the neutrons have a large contribution to the RPA wave 
function of this excitation, the neutron spin current is 
very small, as seen from the short-dashed curves in Fig. 



IV. INTERPLAY OF E2 AND Ml (e,e') CROSS 
SECTIONS 

We use the QRPA transition densities (^5|) for a se- 
lected intrinsic excitation 1+ to calculate its (e, e') cross 
sections with the DWBA code of J. Heisenberg Q]. The 
Ml excitation is fully specified by the transverse tran- 
sition density Jii{r), while the E2 excitation has also 
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a longitudinal (Coulomb) part, P2{r), besides the two 
transverse components J2i^i{r), L' = 1,3. The DWBA 
code uses only J23 ('') ^ input and takes the other trans- 
verse transition density through the continuity equation 
into account. 

The DWBA (e, e') cross sections of the five 1+ excita- 
tions from Table | are plotted in Figs. |^|| and compared 
to experiment. Each one of these 1+ states corresponds 
to the strongest low-energy orbital Ml transition in the 
respective nucleus. These are all the cases in which ex- 
perimental (e,e') cross sections of single 1+ excitations 
in rare-earth nuclei have been published until now. The 
QRPA transition densities for ^^^Gd and ^®^Er, used in 
the DWBA calculations, are shown in Figs. |^ and ^. 

The DWBA cross sections in Figs. ||-|| correspond to 
the scattering angle 9 = 165°, at which the experimental 
data js^,^,^ (dots with error bars) were taken. The 
wave of the incoming electron is distorted by the ground 
state charge distribution of the target nucleus, described 
by a two-parameter Fermi formula. We use for ^^*Sm and 
^^^Gd the experimental values |4^] of the two parameters, 
extracted from electron scattering. Since no experimen- 
tal values are available for the other studied nuclei, the 
values for the latter nucleus are used also for ^^"'■^^^Gd, 
while the experimental values for ^^^Er |Q are employed 
for i68Er. 

It is seen from Figs. H-||that the theoretical Ml cross 
section of the I^K = 1~ 1 excitation (dashed curve) un- 
derestimates systematically the experimental data, espe- 
cially at higher momentum transfer gcff > 0.6 fm~^, i.e., 
q > 0.4 fm~^, corresponding to incident electron energy 
Ei > 40 MeV. This peculiarity was discussed already 
in the Motivation, Sect. I. At small momentum trans- 
fer (close to the photon point) the cross section of the 
accompanying E2 transition to the I^K = 2+1 state 
(dotted curve) is two orders of magnitude smaller than 
the Ml cross section. Thus, the E2 transition will not 
introduce practically any correction to the experimental 
B{M1) value, which is extracted by extrapolating the ex- 
perimental cross section to the photon point q = E^/hc. 

As discussed in Sect. II, the accompanying E2 response 
was identified experimentally only in ^^^Gd as yet [ p^ . 
Its experimental B{E2) value is reproduced in our cal- 
culations together with the corresponding B{M1) value 
(Table I). This is an additional proof for the correct mag- 
nitudes at low q of our theoretical E2 and Ml DWBA 
cross sections for ^'"^^Gd, displayed on the top plot of Fig. 
^. The correct transition probabilities provide, moreover, 
an additional test for our QRPA transition densities P2 (r) 
and J^ii (r) from Fig. |^ because of the simple relationships 



B{E2) 



P2{r)r^dr 





00 



1 r f 1 2 

B{Ml)^-[j^ J,\{rydr . 



(16) 



in magnitude in the region 0.7 < qcs < 0.8 fm ^ 
(0.5 < g < 0.6 fm-i, 50 < Ei < 60 MeV), as seen 
from Figs. |[]|. The sum of the E2 and Ml theoret- 
ical cross sections (continuous curve) agrees very well 
with the available experimental data on (e, e') cross sec- 
tions of single 1+ excitations in rare-earth nuclei. The 
discrepancies in this region, typical for the Ml cross 
section alone, have been removed. It is seen from Figs. 

that the accompanying E2 transition provides an im- 
portant contribution to the measured cross sections for 
0.6 < qes < 0.9 fm-i (0.4 < g < 0.7 fm'^, 40 < E, < 70 
MeV). 

Let us consider further in more detail the 1+ excitation 
in ^^^Gd as an example. Its theoretical Ml cross section 
(dashed line in the top plot of Fig. ^ is plotted once more 
as a continuous curve in Fig. |^ up to a higher momentum 
transfer. It is obtained from the QRPA transition density 
J^ii (r) (continuous line in the corresponding plot of Fig. 
|l|). If one neglects the spin current and uses only the 
convection transition density JTJq (r) (long dashed line 
in Fig. |l]), one obtains the DWBA cross section displayed 
by a dashed curve in Fig. ^. 

The comparison between the continuous and dashed 
curves in Fig. ^ shows that the Ml excitation has a 
predominantly convection character, as it should be ex- 
pected for an orbital transition characterized by the large 
orbit-to-spin ratio Ro.s. in Table |. The total (convection 
plus spin) PWBA cross section (dot-dashed curve in Fig. 
y) is about one order of magnitude smaller than the corre- 
sponding DWBA cross section (continuous curve). Note 
that the PWBA cross section is plotted, as usual, versus 
q, not Soff- 
it is easier to draw qualitative conclusions in PWBA, 
where the cross section can be decomposed into a sum of 
longitudinal and transverse terms (|,^^ : 



(-) 

\dflJ PWBA 

+Vt\\ 



(17) 



Vf = 



cos^(6i/2) 
4 sin-* (61/2)' 



1 + sin2(6'/2) 
8sin''(e'/2) ' 



Vt = Vt/Vt - - + tg2 {9/2) . 



(18) 



The E2 and Ml cross sections become comparable 



Z and e are the nuclear proton number and the unit 
charge, respectively. Note that the squared form factors 
of js],^ have to be multiplied by Att/Z'^ in order to obtain 
our definitions |Q for the form factors |i^L(<z)P (|l7|), co- 
inciding with those of Q . Only the term ^ (q) survives 
in ( p^l ) in the case of Ml transitions, while the E2 exci- 
tations are specified by F^{q) and F^{q). The PWBA 
cross section ( p7| ) is obtained in the relativistic limit (ne- 
glected electron rest mass) and assuming a small excita- 
tion energy in comparison with the momentum transfer, 
i.e., Ej; <C qhc. 
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It follows from that the transverse contribution to 
the E2 cross section is enhanced at 6* = 165° by a kine- 
matical factor (||,|l^) Vr = 58 with respect to the longi- 
tudinal contribution. On the other hand, it is seen from 
Fig. |l| that the amplitude of the longitudinal transition 
density P2if) is one order of magnitude larger than the 
transverse densities J^^j^, (r), L' = 1,3. The squares of the 
corresponding form factors will differ roughly by a factor 
of 100, since the form factors are Fourier transforms of 
the transition densities. Thus, kinematical and dynam- 
ical factors act in opposite directions in backward scat- 
tering and almost compensate each other for 6 = 165°. 

V. QUALITATIVE ANALYSIS OF BACKWARD 
SCATTERING 

More definite conclusions can be drawn from Fig. 
where the PWBA cross section of the I^K = 2+1 excita- 
tion in i56Gd is plotted for 9 = 165°, 175°, 178.5°, and 
179.5°. The total E2 cross section (continuous curve) 
is the sum of the longitudinal C2 (dashed curve) and 
transverse E2 (dotted curve) cross sections, as given by 
(p7|). For small momentum transfer, q « 0.3 fm~^, the 
total E2 cross section at 165° is one order of magnitude 
smaller than the corresponding PWBA Ml cross section, 
shown by a dot-dashed curve in Fig. I Thus, the Ml 
dominance over E2 at small q is preserved in PWBA, 
although it is more pronounced in DWBA where the dif- 
ference is two orders of magnitude at smaller q, as seen 
from the top plot of Fig. ^. 

It is seen from Fig. R that the Coulomb cross section 
C2 (dashed curve) provides the main contribution to the 
E2 cross section at 9 = 165° up to a high momentum 
transfer q < 1.4 fm~^. The transverse E2 cross section 
remains constant in all the four plots because the trans- 
verse factor Vt (|l^) has already reached its maximal value 
at 9 = 165° and does not increase any more appreciably 
with 9. In contrast, the longitudinal C2 cross section 
(dashed curve) is reduced further when going to larger 
backward angles since Vi decreases rapidly. As seen from 
Fig. 1^, the longitudinal and transverse cross sections be- 
come comparable at 175°. 

The total E2 cross section originates exclusively from 
the transverse current at 178.5° in the whole g-range 
studied. As discussed in Sect. II, this is the largest effec- 
tive angle reachable in experiments with backward (e, e') 
scattering Further increase in 9 (179.5°) does not 

produce any change in the E2 cross section. The trans- 
verse diffraction minima occur at higher q-values than 
the longitudinal ones because the transverse transition 
densities i72L'('"), L' = 1,3 are peaked at a smaller radial 
distance than the charge transition density p'^i'r), as seen 
from Fig. |l|. 

The above qualitative analysis in PWBA indicates that 
the important E2 contributions to the experimental Ml 
cross section at intermediate transferred momenta and 



9 = 165°, seen in Figs. ^,|[ originate mainly from the 
Coulomb excitation. Since it is already strongly damped 
at 6* = 178.5°, as seen from Fig. 0, we plot the DWBA 
cross sections for this angle in Fig. |[ It is seen from the 
latter figure that the summed cross section (continuous 
curve) originates mainly from the il/1 excitation (dashed 
curve) for ^cff < 1-2 fm~^. Since both All and E2 cross 
sections are transverse at this angle, the difference in 
magnitude arises mainly from the respective transition 
densities, shown in Fig. |l|. The two cross sections are of 
comparable magnitude for 1.2 < q^g < 2 fm~^. One has 
to expect an E2 dominance for larger momentum transfer 
(not shown in Fig. ^ having in view that lower multipo- 
larities decrease faster for larger q |2j], as discussed in 
Sect. II. 

Therefore, the analysis of the interplay between Ml 
and E2 (e, e') responses of the low-lying orbital 1"'' exci- 
tations leads to the following conclusion. Even at large 
backward angles, where the strong longitudinal C2 re- 
sponse is switched off by the kinematics, the usually ex- 
pected Ml dominance at low q is no more preserved for 
larger momentum transfer, q^g > 1.2 fm"-"^. 

We learned recently that the role of longitudinal 
contributions at 9 = 180° was studied earlier in even- 
odd nuclei |Q. Their importance was demonstrated for 
the ground-state rotational band of ^*^Ta. The multi- 
pole mixing is more important in even-odd nuclei since 
their ground state has a large angular momentum (spin) 
and many multipolarities are allowed by the selection 
rules. In contrast, even-even nuclei have a zero spin 
in the ground state and no multipole mixing occurs in 
transitions to the ground state. Thus, the effect we are 
studying has a different origin: it is entirely due to the 
large moment of inertia of heavy nuclei, as discussed in 
Sect. II. 

The DWBA C2 cross sections were calculated in ^] 
by including the Coulomb and neglecting the transverse 
E2 contributions. The latter, which are strongly en- 
hanced by the transverse kinematical factor at 180°, were 
calculated separately in PWBA for comparison. They 
are smaller than the Coulomb DWBA cross sections 
but, according to the above analysis, they could become 
larger than the longitudinal contribution if calculated in 
DWBA. As noted in |^^, more precise conclusions could 
be drawn after including both longitudinal and trans- 
verse contributions in DWBA calculations, where they 
interfere in contrast to the coherent PWBA sum (fl^) . 
The transitions within a rotational band, studied in p9| , 
are characterized by strong Coulomb excitations because 
the intrinsic wave function does not change in the final 
state and the matrix element of the charge density op- 
erator is large. In our case of intrinsic excitations with 
K"^ = 1+ the Coulomb excitation is due to the predomi- 
nantly orbital character of the low-lying Ml transitions, 
involving a large non spin-flip component in the RPA 
wave function. 

The longitudinal kinematical factor Ve (|l^) vanishes for 
9 = 180°, so that a longitudinal contribution can arise 
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only if the electron rest mass, neglected in ([l7| ), (]i8|), is 
taken into account. The diverging ratio Vr dig ) becomes 
finite in this case liajM : 



qhc 



(A)^ 



(19) 



where mgC^ — 0.511 MeV is the electron rest mass. The 
latter relation in ( pj| ) is valid for <C Ei. It is seen 
from (^9|) that the longitudinal kinematical factor will be 
comparable in magnitude with the transverse one (Vt < 
1) only for unrealistically small momentum transfer q < 
0.005 fm~^ or incident electron energy Ei < 2 MeV. 

A large Coulomb C2 cross section was found however 
in H for 6* = 180° and q^s < 0.7 fm"!. This can be 
understood in PWBA by assuming a transverse electric 
excitation of convection type and using the "virtual pho- 
ton method" j|] for small q to obtain a qualitative esti- 
mate for the relative CL and EL cross sections, given in 
Table XXIV of [i: 



da{CL) 



/'q'hc\'^ 1 /2meC^\' 
v"e7/ \ E^ ) 



(20) 



conv \ Ex ^ Vt ^ Ex 

where also T^t(180°) from (|l|) was used. Note that the 
wrong leading EL term q/M was chosen in not be- 
cause it is spin- flip, as stated in but because it is 
smaller than the correct leading term E^/q for low q. In 
contrast to the severe restrictions (|l^), it is seen from 
(|^) that for any multipolarity L the Coulomb cross sec- 
tion will be larger than the transverse electric one for 
excitation energies E^ < 1 MeV. This is indeed the case 
in 1 49 1 , where the studied first two members of the ground 
state rotational band have very small energies, 0.14 and 
0.3 MeV. 

The situation is more unfavourable in our case, where 
the considered 1+ excitations from Table | have E'^ ~ 3 
MeV. In this case we have still a dominant Coulomb con- 
tribution at 165° for q < 1.2 im-\ but for 9 > 178° 
the C2 cross section is negligible in comparison with the 
transverse E2 cross section, as seen from Fig. |^. This 
result agrees with the qualitative PWBA prediction ( |20| ) 
for a difference of one order of magnitude at 180°. There- 
fore, our results differ qualitatively from those for the ro- 
tational excitations of even-odd nuclei , where the at- 
tention was focused on a fully backward angle, 9 = 180°, 
so that the main effect found there (large C2 contribu- 
tions) is due to the small but non-negligible electron rest 
mass (p^, (pO[). This is not the case with the 1+ exci- 
tations in even-even nuclei, studied in the present work, 
where the Coulomb contribution becomes negligible al- 
ready before reaching a full backward angle. 

Thus, in our case of low-lying orbital 1+ excitations 
the interplay of transverse electric and magnetic contri- 
butions is of prime importance for 9 > 178°. We borrow 
the PWBA estimates in the low-g limit from Tabic XXIV 
of i: 



da{M, L-1) 
da{EL)conv 



e: 



(21) 



where M is the nucleon mass. A nuclear radius i? = 6.5 
fm was assumed for the rare-earth region to obtain the 
latter relationship in (pT|), where the excitation energy E^ 
must be supplied in MeV. The spin contribution to the 
transverse electric cross section is negligible with respect 
to the convection contribution at low q and the ratio 
between the M 1 and E2 cross sections is determined by 
( pT| ) . This relationship should be valid for any scattering 
angle. It applies to both the convection and spin parts 
of the magnetic excitation (M, L — 1). 

It is seen from ( ^ ) that a dominant Ml excitation 
is predicted in PWBA for small rotational energies, e.g. 
Ex < 0.3 MeV but not in our case where Ex 3 
MeV (Table |). The PWBA expectation for a dominant 
E2 transition, which should be valid at least for low q, is 
not supported by our microscopic results, shown in Fig. 
P the Ml cross section (dashed curve) is one order of 
magnitude larger than the E2 response (dotted curve) 
for qcs < 1.2 fm~^. 

The main features of the competing E2 and Ml re- 
sponses to inelastic electron scattering, found by the 
qualitative analysis in this section, are common for all 
the low-lying 1+ excitations from Table |, studied in this 
work. Each of them is the strongest orbital Ml excita- 
tion in the respective nucleus and has the character of a 
weakly collective scissors mode 11,13. 



VI. SUMMARY AND CONCLUSIONS 

We have studied the competition between electric and 
magnetic excitations with K'^ = 1"*" in response to in- 
elastic electron scattering from heavy deformed nuclei. 
This study, whose initial results were reported briefly 
in ||2^,^, is, to our knowledge, the first theoretical in- 
vestigation of possible electric contributions to Ml cross 
sections of even-even nuclei, measured experimentally at 
backward scattering angle. The topic is of some gen- 
eral interest, because the backward electron stattering 
has established itself since three decades as one of the 
main tools for the experimental study of nuclear mag- 
netic dipole excitations due to its particular sensitivity 
towards the magnetic response at backward angles 

Electric contributions have been usually neglected at 
backward angles, apart from a theoretical indication ^] 
for a possible Coulomb excitation in the ground state ro- 
tational band of heavy even-odd nuclei on the example of 
^*^Ta. The multipole mixing is allowed in even-odd nu- 
clei since they have a non-vanishing angular momentum 
in the ground state. In contrast, the mixing is forbidden 
by selection rules for transitions to the ground state of 
even-even nuclei, which has a zero spin. 

However, a rotational band develops in deformed nuclei 
on the top of each intrinsic excitation with K'^ — 1+. The 
Ml transition with I^K = 1+1 is accompanied by an E2 
transition to the first rotational member with I'^K = 2+1 
of this band. The energy separation between the two 
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excitations is small in heavy nuclei with a large moment 
of inertia and comparable with the experimental energy 
resolution in electron scattering. Such a mechanism for 
multipolc mixing was not considered previously, although 
it gains in importance through the systematic occurrence 
of accompanying E2 transitions. It stands in contrast to 
the usual low probability for the random presence of some 
other intrinsic E2 excitation in the vicinity of the studied 
Ml transition. 

The accompanying E2 transition was identified experi- 
mentally only in ^^^Gd as yet ||l^]. We were motivated by 
this work in the search for a possible explanation for the 
systematic deviations of the theoretical Ml form factors 
from experiment at higher momentum transfer, occurring 
only in heavy nuclei, but absent in light ones jllj even 
up to higher transferred momentum. 

We describe the 1^ excitations in heavy nuclei within 
QRPA using a deformed Woods-Saxon potential, BCS 
pairing, and separable residual interactions. Our formal- 
ism was able to provide a good description of 1+ ex- 
citations (see the review Q]) in Titanium ||ll|,^, rare- 
earth | p^ , p^ , ^ and actinidc |l^ nuclei, studied in (e, e'), 
(7,7'), and experiments, e.g., ]T^ , ^ , ^8| -[4^] , re- 

viewed in 0,|l|. 

Using this formalism, we have studied here the com- 
peting E2 and M 1 excitations in all the rare-earth nu- 
clei with experimentally measured (e, e') cross sections. 
These are the strongest Ml excitations in each nucleus. 
They have an orbital nature and correspond to a weakly 
collective low-energy scissors mode @0. The theoret- 
ical energies, B{M1) values, and orbit-to-spin ratios are 
in agreement with experiment. The only available exper- 
imental B{E2) value ||l^ is also reproduced. The nec- 
essary transition densities are obtained microscopically 
from the QRPA wave functions within our formalism, de- 
veloped in a cylindrical basis. The (e, e') cross sections 
are calculated in DWBA from the transition densities. 

The comparison with the experimental (e, e') cross sec- 
tions, measured at scattering angle 9 = 165°, shows 
that the previous systematic deviations of the theoretical 
Ml cross sections from experiment at higher momentum 
transfer q are removed by taking the E2 response into 
account. At low q the theoretical E2 cross section is 
two orders of magnitude smaller than the corresponding 
Ml cross section. Therefore, the E2 response does not 
introduce any corrections to the B{M1) values derived 
by extrapolating the experimental cross sections to the 
photon point (minimal q value). 

The E2 and Ml cross sections become comparable 
in magnitude in the region 0.7 < gcff < 0.8 fm^^ 
(0.5 < g < 0.6 fm-\ < Ei < 60 MeV). The ac- 
companying E2 response provides an important contri- 
bution to the measured cross sections for 0.6 < ^eff < 0.9 
fm-i (0.4 < g < 0.7 fm-\ m < Ei < 70 MeV). The sum 
of the E2 and Ml theoretical cross sections agrees quite 
well with the available experimental data. The qualita- 
tive behaviour, the improvement, and the achieved agree- 
ment with experiment are manifested in a similar way in 



all the excitations studied and exhibit a well pronounced 
systematic character. This allows us to conclude that the 
E2 response of low-lying orbital 1+ excitations to inelas- 
tic electron scattering plays an important role in heavy 
deformed nuclei for higher momentum transfer even at 
backward scattering angles. 

The PWBA analysis of the cross sections has shown 
that the E2 excitation is dominated by the Coulomb C2 
response at 6* = 165° and this qualitative feature is pre- 
served for 9 < 175°. In contrast, the longitudinal C2 
response is already negligible for the largest experimen- 
tally reachable effective backward angle 9 = 178.5° [ p3| . 
It is simply damped by the small longitudinal kinematic 
factor. The interplay between magnetic and electric re- 
sponses is determined at this angle (and beyond it until 
180°) by the competing transverse magnetic and electric 
responses. Both of them receive the main contribution 
from the proton convection current. 

The transverse Ml response is still dominant at 178.5° 
for qcB < 1.2 fm^^, corresponding to incident electron en- 
ergy Ei < 100 MeV. The two responses are comparable 
beyond this energy and the E2 cross section is expected 
to be even the dominant one for Ei > 200 MeV. There- 
fore, the commonly expected Ml dominance in backward 
scattering does not take place at high momentum trans- 
fer even when the full backward angle is reached. This 
is true at least for the orbital Ml excitations studied in 
the present work. 
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APPENDIX A: CYLINDRICAL BASIS 

A large number of formulas for multipole transition 
densities has been derived in the past using non-trivial 
relationships from the differential tensor calculus, see, 
e.g., 1^ and references therein. They are very useful for 
microscopic calculations, but almost all of them are for- 
mulated in a spherical basis, whose properties are essen- 
tial for obtaining the final results. The spherical basis 
is not suitable for deformed nuclei. We follow the ap- 
poach of our former work p3| and derive the necessary 
expressions for deformed basis in cylindrical coordinates 
r = (p, z,ip), where p and if are the azymuth radius and 
angle, and z is the symmetry axis of the nucleus. 

We are using a deformed, axially symmetric Woods- 
Saxon potential pl| ], whose eigenfunctions are expanded 
in the basis of the eigenfunctions of the deformed, axially 
symmetric harmonic oscillator 
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= ^(np, |A|, p) z) ^(A, if) x(S), (Al) 

where j = {np,nz, A,Yi) denotes a set of four quantum 
numbers fully specifying the wave function. x(S) is a 
two-component spinor with projection S = ±1/2 of the 
spin S on the quantization z-axis. The principal quantum 
number Af, parity tt, projection A of the orbital angular 
momentum L and projection K of the total angular mo- 
mentum J = L -|- S on the nuclear symmetry axis z are 
related in the following way: 



TV = 
fs: = A + E, 



2n„ 



|A|, 



7r = (-l> 



(A2) 



The remaining functions in (Al) are expressed ]5l| ] 
through the Hermite and associate Laguerre polynomi- 
als H„, , lLJ : 

^K, |A|,p) = ^fip) = Nll^lV2R^,?^exp(-|)Llfl(,y), 

e 



(A3) 



where shorter notations are introduced for the wave func- 
tions with j denoting either Up or , to remind that both 
of them belong to the same set of four quantum numbers 
specifying the eigenfunction (Al). The norms N and the 
coordinate scaling factors R are given by 

nL'^J = ^np\/{np + \A\y., N„^ = (^/^2"^ njyK 



R± = y/Muj±/h, = /MwIT^- (A4) 

M is the nuclcon mass and lj±_ , ojz are the two frequences 
of the axially symmetric harmonic oscillator. The time- 
reversed basis wave functions ( |Al|) have the form 

*j(r) =f*,(r) - (-l)i/2-^*(np,n,,-A,-S,r). 

(A5) 

Let us introduce short notations D for the derivatives 



of (A3) with respect to p and z: 



d_ 



i^{np,\A\,p)^Df{p) 



|A| 



-0(np, |A|,p), 



Dfip) = -Rj 



;^(np-i,|A| + i,p) 



|A| + lV(np,|A| + l,p) 



d 



dz 
R^ 

d 

— V'(A, ip)^iAtl^{A, If) 



nz ip{nz -1,2;)- Vnz + lip{nz + l,z) , 

(A6) 



The eigenf unctions $i of the axially symmetric Woods- 
Saxon potential ||3l|] are linear combinations of ( Al ) with 
definite projection K and parity vr. Different functions 
are labelled with i, because there are many orthogonal 
functions with the same if. For K = Ki> 1/2: 

$,(i^^r) =;^q^',(np,n,,A,S = i,if,r) 

j 

+ E Cf^An'p, <, A' = A + 1, E' = -i, X, r), 
i' 

A = i^-l/2, (A7) 

where Cj, j = (np^riz), are expansion coefficients found 
after numerical diagonalization of the hamiltonian in 
som e truncated basis with M < A/'max.- The expansion 
( |A7| ) is separated explicitly in two terms with constant 
(A, E), because this property will be used below in the 
analytical derivation of expressions for transition densi- 
ties. The time-reversed state corresponding to (A7) is 
easily obtained from (A5): 



$,(/f , r) = ^ q^jirip, nz, -A, E = -i, -K, r) 

j 

- J2 C]> {n'p. < , -A - 1, E' = i, ~K, r). (A8) 



We describe intrinsic excitations with /^K = \Kf — 
Ki\ = 1. Thus, time-reversed states appear only in 
matrix elements between initial and final states with 
Ki = -1/2, Kf = 1/2 and A = 0. It is explained in 
Sect . D h ow they can be obtained from the expressions 
(^), (D8) by some phase changes. 



APPENDIX B: SYMMETRIES OF OPERATORS 

The charge and current density operators (0J|) have 
the following properties under time-reversal T and her- 
mitian conjugation: 

p\v)=p{v), fp{r)f-' = p{r), 
jt(r)=J(r), Jt(r) = (-l)^J_p(r), 
f J(r)f-i = -J(r), f Jp(r)f-i = (-1)^+^ J_p(r), 



Jo — Jz-, Ju 



-^[X+My), M = ±l, (Bl) 



where J^, = 0,±1 are the cyclic (te nsorial) compo- 
nents of the vector. It follows from ( |Bl[ ) that the multi- 
pole density operators (^ have the properties 

pLmM = i-irPL,-Mir), 
fpLM{r)f-^ = {-lY' PL,-M{r), 
jluMir) = {-lf+'^'+''+'JLL'.-M{r), 
fjLL'M{r)f-' = (-l)^+^'+*^iLL',-M(r). (B2) 
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The m-symmetrized operators change only a 

phase under hermitian and time conjugation. Let us de- 
note with c^, c""" the corresponding phases of such an 
operator B{m): 

B\m) = c"B(m), fB{m)f-^ = c^B{m). (B3) 



It foUows from (B2) and corresponding relationships for 
Q{m)^ J(m) that the m-symmetrized operators ( |l^ ) 
have the properties hsted in Table ||. 

The spin and orbital angular momentum operators are 
symmetrized in the same way as J{m) in (|l^. It is 
seen from Table || that all transition operators have real 
matrix elements in our cylindrical basis ( Al ) , ( A7 ) , ( [A^ ) , 
apart from the current density operators JTll' ('ti, r), 
whose matrix elements are imaginary. Operators with 
771 = +1( — 1) arc hermitian (antihermitian), apart from 
the angular momentum and spin operators in the last 
column of Table ||, which have the opposite property. 

All the operators acquire a phase of the type c^c""" = 
±777. This phase can always be chosen to be -I-777 by 
redefining B{m) B{—m), if necessary. The convention 



C^cT 



(B4) 



has already been taken in the definitions (|4|) into ac- 
count, as seen from Table ||. It can be shown that, after 
adhering to ( B4 ) and neglecting in QRPA the scattering 
terms a^a, ocr, any operator with the properties (B3) 
has the following quasiboson representation: 

t.f>0 

X [A^if, m) + (7 mA{if, m)] 
+b„ ifi, m) [At {if, m) + a mA{if, m)] } , (B5) 



Basis operators (B8) with different tti- values commute 
and determine the same property of the quasiboson op- 
erators (^), 

[A(7/,77l), A\if,m')]=Q, 
[B„{m), Bl{m)]=0, 
[Brj'{m'), Bcr{m)] = 0, for 777' ^ m, 

[B,{m), B,,{m)]^Q, ior a' ^ a. (B9) 

Creation and annihilation operators with the same a, m 
commute as well, bec aus e the hermitian conjugation pro- 
duces only a phase (B3). This property is a direct con- 
sequence from the fact that the QRPA vacuum r^(m)|) 
(|Cl[ ) is not a vacuum for the quasiboson operators ( |B5| ), 
i.e., B^{m)\) ^ 0. Only operators with the same 777 
and different a (longitudinal and transverse) have a non- 
vanishing commutator, reflected in the last relation (|BS|). 

It is seen from Table || that the quadrupole and charge 
density operators are characterized by cr = 1 (longitudi- 
nal type), while the spin, angular momentum, and cur- 
rent density operators, corresponding to axial vectors in 
the (x, ?/)-plane, have a — —1 (transverse type). This 
fundamental property can not be changed by redefining 
the 777-symmetrization in (|lj). The q.p. matrix elements 
of Q{m) and J(rr7) are given in ||l^. Note that we are 
using here a different definition of 7(777) ji^Js^. The q.p. 
matrix elements (B6) of the density operators ( p^ ) have 
the form: 

P2{fi,m,r) = --^U{f,i){f\p2{m,r)\i), 

P2{fi,m,r) = ■^U{f,i){f\p2{m,r)\i), 
U{f,i) = UfVi + u^Vf, V{J,i) = UfVi - UiVf, 



b„{Ji,m) = ~-j={ufVi + aUiVf){f\B^{m)\i), 

b„{fi,m) = -^{ufVi + (TUiVf){f\Bcr{m)\i), (B6) 

(B7) 



a - c™ ^ c^c" - ±1, 



where it, v are the BCS occupation probabilities and 
the indices {i, /) denote the Woods-Saxon eigenfunctions 
(A7),(|A^), used to calculate the single-particle matrix 
elements in (B5). The summation runs only over states 
with K > because time-reversed states (with K < 0) 
are treated explicitly in (B5). 

The q.p. basis is provided by the 7r7-symmetrized 2q.p. 
creation and annihilation operators |lq| . appearing in 

(Pi, 



1 



(i f , 777) — — ;= (a\at + ma 



]4) 



(B8) 



jLL'{hm,r) = l=Vif,i){f\jLL'{m,r)\t). (BIO) 



The q.p. matrix elements ( BIC ) have the same sym- 
metries as those of the remaining operators (pl|). Thus, 
one can write in general for the q.p. matrix elements of 
any operator Ba-{fn) (p^), among those listed in Table ||, 

B{fi, 777) = mB{if, m) — mB{fi, 777) — B{if, 777), 
B{fi,m) — mB{if,m) 

= ~mB{fi, 777) = ~B{if, m), (Bll) 

where the phase a and the radial dependence of density 
operators are skipped for simplicity. The definitions of 
the 777-symmetrized operators ( [T^ have been specially 
chosen to ensure the above uniform properties of their 
matrix elements, valid for both longitudinal (u = 1) and 
transverse (tr — —1) operators. 
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APPENDIX C: MULTIPOLE TRANSITIONS 

The QRPA phonon cre ation operators are ex- 
panded over the q.p. basis 



t.f>0 

-(l>u{fi,m)A{iJ,m) 
+XMh ni)A\if, m) - ti^ifi, m)A{if, m)] , (CI) 

where ip, X are forward- and cf), /i backward-going RPA 
ampHtudes, found numericaUy for each excited state 
h' by solving the RPA equations of motion. The m- 
symmetrized wave function is obtained in the laboratory 
frame after projecting the intrinsic excitation rj,(m)|) on 
a good angular momentum One obtains with this 

wave function the reduced matrix element of any transi- 
tion operator (irreducible tensor of rang L) between 
the nuclear ground state and the considered ^-th RPA 
excitation with a total angular momentum /, its projec- 
tion K on the nuclear symmetry axis, and parity tt: 

{rK = L+l,iy\\ fl II 0,g.s.) 

= ^ E [(ir.Mrlil) 



m=±l 

+ (-i)^+V(|r.(m)fl^_j)' 



(C2) 



The multipolarity L of in ( p2| ) specifies the corre- 
sponding m-symmetrized operators from (|l^), i.e., L = 1 
refers to J^ii{m,r), J{m), L{fn)^ and S{m), while L = 2 
means /52(to, r), J2L' (w, r), or Q{m). The necessary 
cyclic components appearing in (|C^), are easily 

obtained from (p^, where all the operators are taken in 
the quasiboson representation (B5): 



tIm = Ba{m = -cr) + M B^{m = a), B ^ J, 



for B(m) = Jim). (C3) 



The matrix elements in ( |C2[ ), calculated with the help 
of (|C3|), have the form: 

{T^{m) B^{m,r)) ^ ([r^(TO), B^{m,r)]) 
[Kifh'm)baifi,m,r) 



ij>0' 



(C4) 



where b„{fi^m^r) are the q.p. matrix elements of 
Ba{m,r) (B6). They are given by ( B1C| ) in the case of 
transition density operators. The coefhcients F^j fi , m) 
are linear combinations of the RPA amplitudes (CI ): 



Fl^{fi,m) = \y{fi,m) + am ^iy{fi,m). (C5) 



A different notation, F^{fi,m) = F^^{fi,m), G^{fi,m) 
= F'iAfi,m), was used in our previous works, e.g. 
]l6| , ^ 12 1 . It can be shown that the amplitudes ( |C5| ) have 
exactly the same properties as the q.p. matrix elements 
The definitions 



(Bll). The definitions ([IJ) ensure that these properties 
are the same for the q.p. matrix elements of all the op- 
erators ( p^ . One can show i n th is way that the intrinsic 
transition matrix elements (C4) do not depend on the 
signature m, i.e.. 



(r,(m = +l)B,(TO = +l,r)) 
(r^(TO = -l)B^(m = -l,r)). 



(C6) 



Using (C3) and (C6), one obtains for the reduced transi- 
tion matrix elements of the type (C2): 



{rK ^L+l,,y\\fL\\ 0,g.s.: 



0, 



II Tl II 0,g.s.) = 2(r,(TO) i3.(m,r)), (C7) 



where the intrinsic matrix element in the r.h.s of ( |C7| ), 
given by (C4), (C5), can be calculated cither for m = 1 
or for m = —I. It is seen from (C7) that in our m- 
symmetrized formalism the QRPA ground state is not 
only a vacuum for the QRPA annihilation operators ( CI ) , 
r^{m)\) = 0, but also for the transition operators T' 
in the laboratory frame. This property is not trivial, 
because it is not fulfilled for the corresponding transition 
operators Baifn^r), B^{m,r) ( [l^ ) in the intrinsic frame: 
both of them have a non-vanishi ng a ction on the QRPA 
vacuum, as it can be seen from (C4) having in view the 
hermitian properties (p^) , listed in Table ||. 



Half of the terms in microscopic sums of t he ty pe (C4) 
are redundant because of the relationships (p311|) . Thus, 
we make a transition from m-symmetrized (|14[) to the 
usual step-up operators J+, S+ ( [Til) or tensors with 
M = -1-1 for the remaining operators (|14D. Let us define 
with them the following q.p. matrix elements: 



P2l{f^,r)^U{f,^){f\p2l{r)\^) 
P2i{flr)^U{f,i){f\p2i{r)\I) 
JLL'i{ft,r)^V{f,i){f\JLL'i{rm 
JLL'i{flr)^Vif,t){f\JLL'i{r)\i) 



(C8) 



V{f,i){f\J+\t). 
-^V{f,i){f\J+\l) 



The corresponding matrix elements of Q21 are defined in 
analogy with p2i above, and those of and 5+ - in 
analogy with J_|_. The symmetries of such step- up and 
step-down matrix elements can be derived from (Bll) 
using (|l^). These symm etri es allow us to manipulate 
the microscopic sums in ( |C^ ) in order to get rid of the 
redundant terms. Upon insertion in (C7), one arrives to 
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the the following final expression for the r.m.e. ( |C7| ) in 

terms of step-up matrix elements: 

(Ln,z.||fl||0,g.s.) = i Yl [Kif^)bAfhr) 

0<i<f 

+F:^{fi)bMlr)], (C9) 

where 6cr(/z, r), ba{fi, r) are the step- up matrix elements 
( |C8| ) of the multipole transition density operators. The 
notation Q < i < f indicates that the final state / has a 
larger ii"- value than the initial state i. The amplitudes 
of the QRPA phonon wave function are contained in the 
new step- up factors (|C9|), which are related with the 
old ones (05) through: 



- V2[F:ift, -1) - F^ifi, +1) 
F:{fl) = V2[F:{fl +1) - F^ifl -1) 

= -V2[f;(/z, -1) + +1)" , 



(CIO) 



i.e., they are obtained from the old factors (C5) by sim- 
ply replacing q.p. matrix elements of m-symmetrized 
operators with those ( p8| ) of the corresponding step-up 
operators. 



APPENDIX D: MATRIX ELEMENTS OF 
TRANSITION DENSITIES 

We are g oing to derive expressions for the matrix ele- 
ments (|C8 ) of the transition density operators P2i{f) and 
JhL'iir) (13|). They are necessary for the ca lcula tion of 
the E2 and Ml transition densities (15a), (15b). The 
vector scalar product in ( [T^ ) is expressed by the covari- 
ant components of the involved vectors |Q : 

[YLL'Mm\. - (-i)^(^', M + 1, -^i\LM) 

YLL'Aim . Ti = ^(L', M + ^, 1, ^^i\LM) 

YLL.M(f^)»[Vl xTi]i 

-iV2y^{L',M + /i, l,-fi\LM){l,iy, 1, A|l/i) 



;>L'3f+,,(f^)Vi,Ti,A, (Dl) 



where the round brackets denote the Clebsch-Gordan co- 
efficients. 

The action of 0-dependent operators on the basis func- 
tions V'(-Ai 0) (A-1) and spin components on the spinors 
x(S) is calculated analytically. Thus, the integration 
over the azymuth angle (f) and the spin variables, per- 
formed analytically, results in selection rules for the 



quantum numbers A, E of the invol ved Woods-Saxon ini- 
tial and final wave functions (/, ?) (A7). Therefore, only 
the first two basis functions, i.e., the Hermite and La- 
guerre polynomials from (|Al|) , remain in the final expres- 
sions, which have to be integrated over the polar angle 9. 
Due to the different selection rules, we treat explicitly the 
two differe nt t erms with A, A' in the Woods-Saxon eigen- 
functions ( A7), ([A^). Let us introduce a shorter notation 
A^"^ for the product of these two functions, belonging to 
the initial and final Woods-Saxon states (A7), involved 
in the considered matrix element (/,«). The derivatives 
of such an expression with respect to p and z can also be 
written afterwards in a compact form: 

A/'(A/,A„p,z) = ^C{(A/)q(A,) 

kj 

x^pt'ip)i'f'ip)MzHj{^), (D2a) 



Bf\Af,A„p,z)^^Af^ 
opt 



lA,; 



= ^c/(A;)c;(A,)v^(p)z?f (p) M^m^), 



Cf'^{Af,A,,p,z) = ^Af''{Af,A,,p,z) 
= Y,Ci{Af)C]{A,) ^jt'{p)^p{p) M^)D,{z), (D2b) 

kj 



i?/'(A/,A„p,z) = A^/^ ^^c((A/)q(A,) 

^ kj 



iP) Df' ip) + Dt' ip) {p)\ ^k {z)i,, (z) 



Zf\Af,M,p,z) = -^^Af\Af,A,,p,z) 

^Y.^i{Af)C}{A,)i;t'{p)^f'ip) 
kj 

X [Mz)Dj{z) + Dk{z)i^j{z)]. (D2c) 
The factors C), Cl in the above expressions are the ex- 



pansion coefficients (A7) of the initial and final Woods- 
Saxon states, respectively. The symbolic notations d/dpi 
and d/dzi, appearing in (D2b), indicate that only the 
basis function belonging to the initial state is subjected 
to differentiation. It is seen from (D2b), (D2c) that all 
derivatives are expressed through the elementary deriva- 
tives D^{p) and Dj(z) ( |A^ ) with respect to p and z. 
Using the above short notations for derivatives and 



(Dl) for the vecto r sc alar product, one can write the q.p. 
matrix elements (C8) of the multipole transition densi- 
ties (15a), (151) in a more compact form. The matrix 
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elements of the charge transition density (15a) are given 
by: 

P2i{fhr) = eeU{f,i) J Y2i{m\{r)<i>,{r)dn 

X [Af\A + 1, A, p, z) + ^■''^(A + 2, A + 1, p, cos6i, 

p = r sinO, z — r cos9, 



P2i{fi,r) = -esU{f,i) 



X [Af'{l,0,p,z) - Af\0,-l,p,z)]dcoae, 




(D3) 



where $(r) are the Woods-Saxon wave functions (A7). 
The integration over the polar angle in ( p3| ) and similar 
expressions below is limited to the interval < 6 < tt/2 
(or < cos 6* < 1) due to reflection symmetry of the basis 
with respect to the plan e z = 0. 

The matrix elements ( |C8| ) of the transverse (current) 
multipole density operators J^LL'i{r) ( [T^ ) can be decom- 
posed into convection and spin parts, originating from 
the convection and magnetization current density opera- 
torsj^(r) (l)andj^(r) (0), 

JLL'i{n,r) = jEL'i{fhr)+jiLn{fhr), (D4) 

and the same for J7ll'i(/*, All of them are purely 
imaginary in our basis, see Table They are neces- 
sary for the calculat ion of the the transverse (current) 
transition dens ities ( [l5b| ). The matrix elements (D4), 
obtained from (Dl) and expressed in nuclear magnetons 
Pn = eh/{2Mc) w 0.1052 efm, have the form: 
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JL'iU^, r) = -ifiNg'VU, i)V2j2{L', 1 + 1, -fi\Ll){l, V, 1, A|l, -y) f YL',i+^{m,, [<^>\{r)S,,xMr)] dQ. (D6) 

\ 



One obtains from (D5) the convection matrix elements of the transverse All current with LL' =11 and E2 currents 
with LL' =21, 23: 



Jni{fhr) = -iefiNV{f,z)\j-^ ^ 



Bf\A' + l,A',p,z) + 



(|A'|-A- 



Af'{A' + l,A',p,z) 



C^*(A' + 1, A',p,2)|dcos6', (D7a) 



B^'iA' + l,A',p,z) + 



(|A'|-A') 



Af\A' + l,A',p,z) 



+^C^'{A' + l,A',p,z) \dcos0 



(D7b) 



J2'ii(/*,'^) 



lepN 



A'=A,A+1 



(z2 - V) 



,2\ r 



zB-^^(A' + 1, A', p, 2) + pCf\A' + 1, A', p, z) 



zA 



-^■'^VA' + l,A'.p.z)l<icos6l, 
J 



(D7c) 



where A = — 1/2 in (D7) is determined by the i^- number of the initial state Ki. Let us note that here and below 
we have always Kf — Ki + 1. Although A > in the above two-term sums over A', the expressions are written 
in a more general way in order to use them also fo r th e corresponding matrix elements J^LL'i^f^T''')^ involving a 



time-reversed initial state. They are obtained from (D7) by putting A — Ki — 1/2 = —1 and multiplying the first 
term (A' = —1) with an overall factor - 1, i. e., the two-term sums run over A' = —1,0 in this special case. The 
same structure was written explicitly in (D3) for p2i{fi,r). The phase -1 originates from the second term of the 
time-reversed Wood-Saxon wave function (A8). 



The spin matrix elements of the transverse Ml and E2 currents, obtained from (D6), have the form 



r) = -tgS^^vif, [^^'(A + 1, A, p, z) - Rf\A + 2, A + 1, p, z) 

+ - [Af\A + 1, A, p, z) - Af'iA + 2, A + 1, p, z)] - 2Zf\A + 1, A + 1, p, z) 
P 



Rf\A + 2, A, p, z) - Rf'{A + 1, A + 1, p, z) + -A^^A + 2, A, p, z) 

P 



d cos9. 



(D8a) 



1 /3 



Rf\A + l,A,p,z)- Rf\A + 2,A + l,p, z) 



^- [Af\A + 2, A + 1, p, z) - Af\A + 1, A, p, z)] + 2Z-''*(A + 1, A + 1, p, z) 
P 



Rf\A + 1, A + l, p,z)- Rf\A + 2, A, p, z) + -^■''^(A + 2, A, p, z) 

p 



d cos9, 



(D8b) 
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Rf%A + 1, A, p, z) - Rf'{A + 2, A + 1, p, z) 



p[Rf\K + 2, A, p, z) - Rf\K + 1, A + 1, p, z)] - 2A^'(A + 2, A, p, z) 

+ - W\K + 1, A, p, z) - A^'(A + 2, A + 1, p, z)l 
P 

(2z2 -3p2)Z-''*(A + l,A + l,p,;z) -5p2z^'(A + 2,A,p,z)] ldcos6i, 



The two terms with A' = A, A + 1 from (D7) are writ 



ten exphcitly in the expressions (D8) because their sym- 
metry is lost in t he s pin matrix elements. Nevertheless, 
one can still use (D8) also for the corresponding matrix 
elements Jiuiif^^''^)- They are obtained from ( ]D8D by 
putting A = Ki — 1/2 — —1 and multiplying with -1 
each functional A, R, or Z (D2), which has -1 as a s ec- 
ond argument. The radius r in the r.h.s. of (D7), (D8) is 
a constant with respect to the integration. Thus, p and 
z are not independent integration variables but simple 
functions of the integration angle 6, given in (|D3|). 



(D8c) 



17 



* Permanent address: Institute for Nuclear Research 
and Nuclear Energy, Bulgarian Academy of Sci- 
ences, BG-1784, Sofia, Bulgaria. Electronic address: 
roland.nojarov@uni-tuebingen.de 

^ Electronic address: amand.faessler@uni-tuebingen.de 

* Electronic address: michael.dingfelder@uni-tuebingen.de 
[1] G. A. Peterson and W. C. Barber, Phys. Rev. 128, 812 

(1962). 

[2] W. C. Barber, J. Goldemberg, G. A. Peterson, and Y. 

Torizuka, Nucl. Phys. 41, 461 (1963). 
[3] W. C. Barber, Ann. Rev. Nucl. Sci., 12, 1 (1962). 
[4] L. W. Fagg, Rev. Mod. Phys. 47, 683 (1975). 
[5] H. Uberall, Electron Scattering from Complex Nuclei, 

(Academic Press, New York, London, 1971). 
[6] S. Raman, L. W. Fagg, and R. S. Hicks, Internat. Rev. 

Nucl. Phys. 7, 355 (1991). 
[7] R. Nojarov, Progr. Part. Nucl. Phys. 34, 297 (1995). 
[8] N. Lo ludice and F. Palumbo, Phys. Rev. Lett. 41, 1532 

(1978); Nucl. Phys. A326, 193 (1979); 

G. De Franceschi, N. Lo ludice, and F. Palumbo, Phys. 
Rev. C 29, 1496 (1984). 

[9] N. Lo ludice and A. Richter, Phys. Lett. B228, 291 
(1989); B304, 193 (1993); 

N. Lo ludice, F. Palumbo, A. Richter, and H. J. Wortche, 

Phys. Rev. C 42, 241 (1990). 
[10] A. Richter, Nucl. Phys. A522, 139c (1991); A553, 417c 

(1993); Progr. Part. Nucl. Phys. 34, 261 (1995). 
[11] R. Nojarov, A. Faessler, and P. O. Lipas, Nucl. Phys. 

A533, 381 (1991). 
[12] A. Faessler, R. Nojarov, and F. G. Scholtz, Nucl. Phys. 

A515, 237 (1990). 
[13] R. Nojarov, A. Faessler, P. Sarriguren, E. Moya de 

Guerra, and M. Grigorescu, Nucl. Phys. A563, 349 

(1993). 

[14] T. Guhr, H. Diesener, A. Richter, C. W. de Jager, H. de 
Vries, and P. K. A. de Witt Huberts, Z. Phys. A 336, 
159 (1990). 

[15] A. Richter, Nucl. Phys. A507, 99c (1990). 

[16] R. Nojarov and A. Faessler, Nucl. Phys. A484, 1 (1988). 
[17] R. Nojarov and A. Faessler, Z. Phys. A336, 151 (1990). 
[18] D. Bohle, A. Richter, K. Heyde, P. Van Isacker, J. 

Moreau, and A. Sevrin, Phys. Rev. Lett. 55, 1661 (1985). 
[19] S. Raman, C. H. Malarkey, W. T. Milner, C. W. Nestor 

Jr., and P. H. Stelson, At. Data Nucl. Data Tables 36, 1 

(1987). 

[20] R. Nojarov and E. Nadjakov, Nucl. Phys. A397, 29 

(1983). 

[21] D. Bohle, A. Richter, U. E. P. Berg, J. Drexler, R. D. 
Heil, U. Kneissl, H. Mctzgcr, R. Stock, B. Fischer, 

H. Hollick, and D. Kollewe, Nucl. Phys. A458, 205 
(1986). 

[22] S. Dobert et aJ., Verhandlungen der DPG 6, 747 (1995). 
[23] M. Petraitis, J. P. Connelly, Hall Crannell, L. W. Fagg, 

,J. T. O'Brien, D. I. Sober, .]. R. Deininger, S. E. 

Williamson, R. Lindgren, and S. Raman, Phys. Rev. C 



49, 3000 (1994). 
[24] A. S. Reiner and L. J. Tassie, Nuovo Cim. 37, 1524 
(1965). 

[25] F. W. Hersman, W. Bertozzi, T. N. Buti, J. M. Finn, 
C. E. Hyde-Wright, M. V. Hynes, J. Kelly, M. A. Ko- 
vash, S. Kowalski, J. Lichtenstadt, R. Lourie, B. Mur- 
dock, B. Pugh, F. N. Rad, C. P. Sargent, and J. B. Bel- 
licard, Phys. Rev. C 33, 1905 (1986). 

[26] R. S. Willey, Nucl. Phys. 40, 529 (1963). 

[27] R. Nojarov, A. Faessler, and M. Dingfelder, J. Phys. G 
20, Llll (1994). 

[28] M. Dingfelder, R. Nojarov, and A. Faessler, Progr. Part. 
Nucl. Phys. 34, 319 (1995). 

[29] R. Nojarov, A. Faessler, and M. Dingfelder, Phys. Rev. 
C 51, 2449 (1995). 

[30] A. Faessler and R. Nojarov, Phys. Rev. C 41, 1243 
(1990). 

[31] J. Damgaard, H.-C. Pauli, V. V. Pashkevich, and V. M. 

Strutinski, Nucl. Phys. A135, 432 (1969). 
[32] Y. Tanaka, Y. Oda, F. Petrovich, and R.K. Sheline, Phys. 

Lett. B83, 279 (1979). 
[33] K. A. Erb, J. E. Holden, I. Y. Lee, J. X. Saladin, and 

T. K. Saylor, Phys. Rev. Lett. 29, 1010 (1972); 

H. J. WoUcrsheim and Th. W. Elze, Nucl. Phys. A278, 

87 (1977); 

H. Fischer, D. Kamke, H. J. Kittling, E. Kuhlmann, 
H. Phcht, and R. Schormann, Phys. Rev. C 15, 921 
(1977); 

E. Schormann, H. Fischer, and E. Kuhlmann, Phys. Rev. 
C 16, 2165 (1977). 
[34] S. Raman, C. H. Malarkey, W. T. Milner, C. W. Nestor 
Jr. and P. H. Stelson, At. Data Nucl. Data Tables 36, 1 
(1987). 

[35] A. H. Wapstra and G. Audi, Nucl. Phys. A432, 1 (1985). 
[36] M. I. Baznat and N. I. Pyatov, Yad. Fiz. 21, 708 (1975) 

[Sov. J. Nucl. Phys. 21, 365 (1975)]. 
[37] P. Sarriguren. E. Moya de Guerra, R. Nojarov, and A. 

Faessler, J. Phys. G 19, 291 (1993); 20, 315 (1994). 
[38] D. Bohle, A. Richter, W. Steffcn, A. E. L. Dicperink, 

N. Lo ludice, F. Palumbo, and O. Scholten, Phys. Lett. 

B 137, 27 (1984). 
[39] U. Hartmann, D. Bohle, F. Humbert, and A. Richter, 

Nucl. Phys. A499, 93 (1989). 
[40] C. Djalaii, N. Marty, M. Morlet, A. Willis, J. C. Jour- 

dain, D. Bohle, U. Hartmann, G. Kiichler, A. Richter, 

G. Caskey, G. M. Crawley, and A. Galonsky, Phys. Lett. 

B 164, 269 (1985). 
[41] H. H. Pitz, U. E. P. Berg, R. D. Heil, U. Kneissl, R. Stock, 

C. Wesslborg, and P. von Brentano, Nucl. Phys. A492, 

411 (1989). 

[42] D. Bohle, G. Kiichler, A. Richter, and W. StefFen, Phys. 

Lett. B 148, 260 (1984). 
[43] A. Faessler, R. Nojarov, and T. Taigel, Nucl. Phys. 

A492, 105 (1989). 
[44] D. A. Varshalovich, A. N. Moskalev, and V. K. Kherson- 

skii. Quantum Theory of Angular Momentum, (World 

Scientific, Singapore, 1988). 
[45] R. Nojarov and A. Faessler, Phys. Rev. C 49, 213 (1994). 
[46] J. Heisenberg and H. P. Blok, Ann. Rev. Nucl. Part. Sci. 

33, 569 (1983). 



18 



[47] H. de Vries, C.W. de Jager and C. de Vries, At. Data 

Nucl. Data Tables 36, 495 (1987). 
[48] E. Moya de Guerra (private communication). 
[49] M. Nishimura, E. Moya de Guerra, and D.W.L. Sprung, 

Nucl. Phys. A435, 523 (1985). 
[50] H. Theissen, Springer Tracts in Mod. Phys. 65, 1 (1972). 
[51] L. A. Copley and A. B. Volkof, Nucl. Phys. 84, 417 

(1966). 

[52] R. Nojarov and A. Faessler, Nucl. Phys. A572, 17 (1994). 



FIG. 1. QRPA transition densities (^ for Ml {L = 1) 
and E2 [L — 2) excitations of the K'^ = 1"*" state from Table 
I in ^^®Gd. The total (continuous curve) transverse transition 
density JT^^.' ('") ^ sum of the proton convection (long dashed 
curve) and magnetization (short dashed curve for neutrons 
and dotted curve for protons) transition densities. 

FIG. 2. The same as in Fig. § but for ^**Er. 



TABLE I. Comparison with experiment of QRPA excita- 
tion energies E^, upwards (|) Ml and E2 transition prob- 
abilities and orbit/spin ratios Ro.s. jl^ of the Ml matrix 
elements for the strongest orbital = l""" excitation in each 
rare-earth nucleus with experimentally measured (e, e') cross 
section of this excitation [ p8[|39| , ^ . 



Nucleus 


exp." or 


Ex 


B{M1) 
(pl) 


B{E2y 
(e^fm-*) 






theor. 


(MeV) 




i^^Sm 


exp." 


3.20 


0.80(20) 




> 0.8 




th. 


3.16 


0.87 


28 


7.3 


is^Gd 


exp." 


2.94 


0.85(9) 








th. 


2.94 


1.00 


30 


6.8 


«6Gd 


exp.-*^ 


3.07 


1.30(20) 


40(6) 


> 1.3 




th. 


2.90 


1.24 


42 


7.8 


i5«Gd 


exp.^ 


3.20 


0.77(9) 








th. 


3.48 


0.89 


37 


7.0 


i68Er 


exp.'* 


3.39 


0.90(20) 








th. 


3.35 


0.88 


18 


9.6 



° A reference for the experimental and B{Ml) values 
is given in this column. 

Experimental data are available only for ^^^Gd |l|]. 
^ Experimental data from (e,e') and {p,p') [Q. 

Ref. [E2I. 



" Ref. |3 
f Ref. 
3 Ref. 
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FIG. 3. DWBA (e,e') cross sections for Ml (dashed 
curves) and E2 (dotted curves) excitations, corresponding to 
the QRPA state with K"" = 1+ from Table | in ^^^Sm (top 
plot) and ^^*Gd (bottom plot). Scattering angle 6 — 165°. 
For each nucleus the sum of the Ml and E2 cross sections 
(continuous curve), plotted versus the effective momentum 
transfer qcs jisj , is compared to experimental data (dots with 
error bars) from pij and respectively. 



FIG. 4. The same as in Fig. | but for ^^^^'^^^Gd with ex- 
perimental data from jssj^, respectively. The DWBA cross 
sections for ^^''Gd are obtained with the QRPA transition 
densities from Fig. |l| 



FIG. 5. The same as in Fig. | but for ^'^^Er with experi- 
mental data from [^ . The DWBA cross sections are obtained 
with the QRPA transition densities from Fig. |^. 



FIG. 6. The theoretical Ml cross section for ^^®Gd from 
the upper plot of Fig. ^ is plotted here as a continuous curve 
up to a higher momentum transfer geff- Dashed curve: the 
contribut ion fro m th e proton convection transition density 
J^f{r) (p^), (p7aD alone. These DWBA resuhs are com- 



pared also to the total (convection plus spin) Ml cross section 
in PWBA (dot-dashed curve). 



TABLE II. Phases c", c"^ (B3) of the m-symmetrized op- 



erators 
tively. 



(|l4j) under hermitian and time conjugation, respec- 
c is the phase of the matrix elements in the basis 



(A1),(A7),(A8) 



Phase 


Q(m), p2{m,r) 


JLL'{m,r) 


Jim), L{m), S{m) 


c» 


-f m 


+771 


—m 


c-^ 


-\-m 


— m 


+m 


c^ 


+ 1 


-1 


+1 



FIG. 7. Decomposition of the theoretical PWBA electric 
cross section (continuous curve: total E2) into a sum of its 
longitudinal C2 (dashed curve) and transverse E2 (dotted 
curve) components. They correspond to the E2 excitation 
in ^^^Gd, whose DWBA cross section is shown in the upper 
plot of Fig. ^. The PWBA cross section is plotted here for 
scattering angles = 165°, 175°, 178.5°, and 179.5° up to a 
higher momentum transfer Qeff- 



FIG. 8. The same as in the upper plot of Fig. 
for a scattering angle 6 = 178.5°. 



'Gd, but 
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